Chemical Equilibration in Pb+Pb collisions at the SPS 
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Abstract 



An improved statistical model with excluded volume corrections and res- 
onance decays is introduced and applied to the complete presently available 



. set of particle ratios as measured by the various experiments at the SPS in 

m , 

Pb+Pb collisions. The results imply that a high degree of hadrochemical 

ON 

equilibration is reached at chemical freeze-out in Pb+Pb collisions. 
Heavy ion collisions at ultra-relativistic energies are studied to look for signs of the 
production of a quark-gluon plasma phase which subsequently hadronizes. In this context 
one of the crucial questions is whether thermal and chemical equilibrium is achieved at some 
stage of the collision. Applying a statistical model which assumes equilibrium, and testing 
experimental data against model predictions is one way of testing reality against a thermally 
and chemically equilibrated fireball at the point of hadro-chemical freeze-out. 

The present statistical model - like its predecessor which was presented in - is based 
on the use of a grand canonical ensemble to describe the partition function and hence the 
density of the particles of species i in an equilibrated fireball: 

9i f°° P 2 dp m 
U% 2vr 2 io e (^(p)-«)/r ± 1 [ ) 

with particle density rij, spin degeneracy gi, h — c = 1, momentum p, total energy E and 
chemical potential /i, = jisBi — UsSi — Hi 3 If- The quantities B i} Si and If are the baryon, 



strangeness and three-component of the isospin quantum numbers of the particle of species i. 
The temperature T and the baryochemical potential /xg are the two independent parameters 
of the model, while the volume of the fireball V, the strangeness chemical potential /is, and 
the isospin chemical potential /ij 3 are fixed by the three conservation laws [] for 

baryon number: V n^Bi 

i 

strangeness: V riiSi 

i 

and charge: V ritlf 

i 

Here, Z and N are the proton and neutron numbers of the colliding nuclei. The hadronic mass 
spectrum used in the calculations extends over all mesons with masses below 1.5 GeV and 
baryons with masses below 2 GeV. This limits the temperature up to which thermal model 
calculations are trustworthy to T max < 185 MeV. We note, however, that calculations with 
higher temperatures should anyway be considered with caution as the mass spectrum for 
heavier hadrons is not sufficiently well known. Maybe somewhat unexpectedly the neutron 
excess in Pb plays only a minor role in the determination of particle ratios. 

To take into account a more realistic equation of state we incorporate the repulsive 
interaction at short distances between hadrons by means of an excluded volume correction. 
A number of different corrections have been discussed in the literature. Here we choose that 
proposed in |3|,|J: 

p exd iT, = p id -° as (T, p); with p. = n- v eigen p exd iT, /x). (5) 

This thermodynamically consistent approach to simulate interactions between particles 
by assigning an eigenvolume v e i gen to all particles, modifies the pressure p within the fireball. 
Equation [5] is recursive, as it uses the modified chemical potential p, to calculate the pressure, 

1 These conservation laws apply strictly only for quantities which are evaluated over the complete 
phase space. 
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while this pressure is also used in the modified chemical potential, and the final value is found 
by iteration. Particle densities are calculated by substituting \i in eq. [1] by the modified 
chemical potential fi. 

The eigenvolume has to be chosen appropriately to simulate the repulsive interactions 
between hadrons, and we have investigated the consequences for a wide range of parameters 
for this eigenvolume in ||. Note that the eigenvolume is v e i gen = 4^nR 3 for a hadron 
with radius R. Assigning the same eigenvolume to all particles can reduce particle densities 
drastically but hardly influences particle ratios. Ratios may differ strongly, however, if 
different values for the eigenvolume are used for different particle species. 

Two different scenarios were explored: First, we chose the radius of all baryons according 
to the charge radius of the proton, which lies at 0.8 fm, and assigned a smaller radius of 
0.62 fm to all mesons, as suggested in ||. This drastic correction reduces the thermally 
produced particle density in Pb+Pb collisions by a factor of seven as compared to the ideal 
gas case. Furthermore, because baryons take up more space than mesons, their creation 
is suppressed in favor of meson production. Hence the meson to baryon ratio increases 
strongly. This is illustrated in Figure [I| where we plot, for different chemical potentials, the 
temperature dependence of the pion/nucleon ratio and compare it to predictions from the 
ideal gas scenario. 

However, the nucleon-nucleon or pion-pion interaction is not repulsive at such large dis- 
tances. A more physical approach is to determine, for nucleons, the eigenvolume according to 
the hard-core volume known from nucleon-nucleon scattering ||. Consequently, we assigned 
0.3 fm as radius for all baryons. For mesons we expect the eigenvolume not to exceed that 
of baryons. Therefore, to illustrate the effect, we kept the ratio of meson to baryon radii as 
above, implying a meson radius of 0.25 fm. As can be seen in Figure [I], the resulting particle 
ratios are much closer to predictions using an ideal gas scenario; absolute yields are reduced 
by about 30 %. In the absence of detailed information about the meson-meson interaction 
at short distances we assumed for the following calculations that Rbaryon = Rmeson=0.3 fm. 

After thermal "production" , resonances and heavier particles are allowed to decay, there- 
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fore contributing to the final particle yield of lighter mesons and baryons. Decay cascades, 
where particles decay in several steps, are also included. A systematic parameter regulates 
the amount of decay products resulting from weak decays. This allows to simulate the 
different reconstruction efficiencies for particles from weak decays in different experiments. 

This model is now applied to Pb+Pb collisions at maximum SPS energy. We have used 
all data currently available 0. We adjust the free parameters T and \xb such that they 
reproduce best all particle ratios available at the moment. We did not include in the fit the 
20/(7T + + 7r~) ratio, because, for this ratio, the currently available two experimental values 
exhibit a rather large discrepancy. The results for the best \ 2 ( see below) are shown in Table 
|. For technical reasons only one of the p/p ratios (the NA49 value) was included in the fit. 
All experimental particle ratios are taken from data integrated over transverse momentum 
and integrated over rapidity y to the extent data are available as shown in Table [I]. Using 
this procedure strongly reduces the possible influence on particle ratios of dynamical effects 
such as hydrodynamic flow or particle production from a superposition of fireballs 0. 

The criterium for the best fit was either a minimum in 



or a minimum in the quadratic deviation 



2 = v (nT p - - n? odc1 ) 2 



In the above equations 7£ model and 7^ xp ' are the ith particle ratio as calculated from 
our model or measured in the experiment, and Oi represent the errors in the experimental 
data points as quoted in the experimental publications. We have used both the x 2 an d the 
quadratic deviation measure to estimate the influence of possible systematic errors which 



2 The A/A ratio from the NA49 collaboration will be revised (P. Seyboth, NA49 collaboration, 
private communication) and is therefore not included in Table |. For the same reason we have 
replaced the Z~ /A ratio from NA49 by the ratio (H + + E~)/(A + A). 
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are generally not included in the data. The deviation between these two analyses gives an 
indication of the accuracy of the parameters extracted from this data set. 

As one can see from Figure |2], the best fit in terms of \ 2 was achieved at T = 168 ± 
2.4 MeV, \i B = 266 ± 5 MeV, with /i s = 71.1 MeV and /i l3 = -5.0 MeV. The minimal 
quadratic deviation is found at T = 164 MeV, /!# = 274 MeV. These small differences give 
an indication of the systematic uncertainties of the procedure. 

The overall agreement between model and data is quite good, as can be seen in Figure ||| 
and Table |. Choosing a slightly different excluded volume correction with Rb a r y on=0.3 fm 
and R me son=0.25 fm yields very similar results. Using significantly larger eigenvolumes leads 
to much poorer agreement. Comparison between data and the model, e.g., with Rbar y on=0.8 
fm and R meson =0.62 fm yields xLn = ^"i which is roughly 5 times as large as the value 
shown in Figure 0. In any case, we have excluded such large radii for independent physics 
reasons as discussed above. 

Furthermore, using such large radii leads to an eigenvolume of all particles which would 
occupy 55% of the total volume and could therefore not be considered a "correction" . The 
total fireball volume would increase to roughly 20000 fm 3 , exceeding even the fireball volume 
estimated using pion interferometry. As discussed in PDipTJ, the total fireball volume in 



central Pb+Pb collisions at thermal freeze-out should be about 13500 fm 3 . 

In the small excluded volume scenario with Rbaryon=Rmeson=0.3 fm the fireball volume 
of 2800 fm 3 is considerably smaller. The corresponding pion density is then 0.60 pions/fm 3 , 
significantly exceeding the measured pion density of roughly 0.12 pions/fm 3 [[21]]. This is 
not surprising, however, as the calculated pion density of 0.6/fm 3 is determined at chemical 
freeze-out corresponding to T=168 MeV. If one lets this fireball expand isentropically to 125 
MeV, the temperature roughly corresponding to thermal freeze-out as indicated by particle 
spectra |22fl and two-pion correlations [EH] the corresponding pion density is 0.084/fm 3 , close 



to the experimental value. Note that the calculated pion density would increase further by 
about 30 % if one were to reduce to zero the meson radius in the excluded volume correction. 
The present results, in particular those involving multi-strange baryons, imply that no 



separate strangeness suppression factor is needed to describe the available Pb+Pb data at 
SPS energy. In fact, the mean value of the experimental to calculated yields ratios involving 
AS = 1, i.e. those which are sensitive to a possible overall strangeness suppression, is 0.96 ± 
0.05, consistent with unity. This conclusion differs from that reached in a recent investigation 
2B| where, however, only a very resticted set of ratios was used for comparison with thermal 



model predictions. An interesting anomaly would arise if the 0-meson yield converges to 
the low value reported by the NA49 collaboration (see Table |), since this meson carries two 
units of hidden strangeness. To reconcile this with the results by the WA97 collaboration 
on cascade or omega-baryon production would be a challenge. 

We further note that the improved model discussed here was also applied to the AGS 
data collected in [0J. The best fit, obtained for Rbar y on=Pmeson=0.3 fm, yields T = 125 
(+3 - 6) MeV and //# = 540 ±7 MeV, well in line with the calculations reported in [|TJ. 
Here, the corresponding 7r + and proton densities of 0.051/fm 3 and 0.053/fm 3 agree well with 
those estimated from particle interferometry p4] , p5| (0.058/fm 3 and 0.063/fm 3 , respectively) 
implying that, at AGS energy, thermal and chemical freeze-out take place at nearly identical 
temperatures. 

The good agreement between the predictions of the thermal model and the measured 
particle ratios implies that thermal and chemical equilibrium is established (or at least closely 
approached) in the fireball at hadrochemical freeze-out. Furthermore, it is interesting to 
note that the resulting temperature and chemical potential values are very close to where 



we believe is the phase boundary between hadronic matter and the quark-gluon plasma [26 



It is, therefore, quite probable that the system crosses this phase boundary shortly before it 
freezes out hadrochemically. 
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TABLES 

TABLE I. Experimental particle ratios compared to 

model predictions with Rbaryon = R-mcson = 0.3 fm, T=168 MeV, ^£=266 MeV, /Us=71.1 MeV, 
H h =-5.0 MeV. 

( a ) : feeding from weak decays excluded, 
( 6 ) : feeding from weak decays included, 
( c ) : cuts exclude feeding of A from S ± and S . 

In all three cases feeding in the model was tuned accordingly. In all other cases feeding from weak 
decays is assumed to be 50%. 
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model 



exp. data 



cxp. 



y-range 



ref. 



(p-p)/h- 
p/p 
P/P 
d/d 

7r~/"7T + 

K°/h" 
A/h~ 
A/K° 
K+/K - 

A/A 
5-/A 
H+/A 

5+ + H-)/(A + A) 
O/H 



0.238 

0.045 

0.060 

1.78 10~ 3 

1.05 

0.087 

0.137 

0.126 

0.096 

0.76 

1.90 

1.90 

0.102 

0.102 

0.185 

0.228 

0.114 

0.228 

0.228 

0.53 

0.154 



0.228(29) ^ 

0.055(10) 

0.085(8) 0>) 

0.94(27) 10~ 3 

1.1(1) 

0.081(13) 

0.125(19) 

0.123(20) 

0.077(11) 

0.63(8) 

1.85(9) 

1.8(1) 

0.131(17) W 

0.110(10) ( c ) 

0.188(39) 

0.206 (40) < c ) 

0.13(3) 

0.232(33) 

0.247(43) 

0.383(81) 

0.219(45) 



NA49 
NA44 
NA49 
NA44 
NA49 
WA98 
NA49 
WA97 
WA97 
WA97 
NA44 
NA49 
WA97 
WA97 
NA49 
WA97 
NA49 
NA49 
WA97 
WA97 
WA97 



0.2-5.8 

2.3-2.9 

2.5-3.3 

midrapidity 

all 

2.3- 2.9 
all 

2.4 - 3.4 
2.4 - 3.4 
2.4 - 3.4 

2.4- 3.5 
all 

2.4 - 3.4 
2.4 - 3.4 
3.1-4.1 
2.4 - 3.4 
3.1-4.1 
3.1-4.1 
2.4 - 3.4 
2.4 - 3.4 
2.4 - 3.4 



m 




m 
m 
m 
a 


m 
m 















20/(7T+ +7T 



19.0 10~ 3 
19.0 10~ 3 



21(6) 10~ 3 
12.2(13) 10~ 3 



NA50 
NA49 



2.9 - 3.9 
all 
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FIG. 1. The influence of different excluded volume corrections for baryons and mesons on the 
7r + /p ratio. 
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FIG. 2. x 2 an d quadratic deviation for the comparison between model particle ratios and data 
explored over a wide range of parameters. The dot represents the parameter set with minimum 
X 2 , the square the set with minimum quadratic deviation. 
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FIG. 3. Comparison between model and experimental particle ratios. For experimental data, 
errors and information about feeding corrections see Table 1 and references there. 
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